Singularity theory study of overdetermination in models for L— H transitions 
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Two dynamical models that have been proposed to describe transitions between low and high 
confinement states (L-H transitions) in confined plasmas are analysed using singularity theory and 
stability theory. It is shown that the stationary-state bifurcation sets have qualitative properties 
identical to standard normal forms for the pitchfork and transcritical bifurcations. The analysis 
yields the codimension of the highest-order singularities, from which we find that the unperturbed 
systems are overdetermined bifurcation problems and derive appropriate universal unfoldings. Ques- 

. . . tions of mutual equivalence and the character of the state transitions are addressed. 
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It is a well-known fact that an overdetermined system of equations usually has no exact solutions. In this paper 
O . we report a novel application of singularity theory methods to resolve a subtle case of overdetermination in two 
O I dynamical systems that model L(low)-II(high) confinement state transitions and associated edge-localized modes 
1^ . (ELMs) in confined plasma devices ||,^ . The analysis also addresses the much-discussed question of whether second- 
VO ' order state, or "phase" , transitions occur in these systems. Since both models are based on sound physics and seek to 
describe the same phenomena, we discuss briefly the issue of equivalence, in terms of the singularity theory results. 

The semiotics and dissemination of singularity theory owe much to the elementary catastrophe theory proposed by 
Thom [Q. In substance, however, the provenance of singularity theory can be traced to the work of Poincare, and the 
original exposition was by Whitney [|. It was subsequently developed rigorously and extended by many others, e.g. 
I0'il0l' Successful applications have included diverse problems in mechanical, biological, and chemical systems. 
OQ ] This is the first systematic application to bifurcation problems in plasma physics. 

In the singularity theory approach, the qualitative properties of a dynamical system are characterized by classifying 
the singularities in the set of stationary states, or bifurcation diagram, over the parameter space. In the bifurcation 
diagram of an idealized dynamical model, a degenerate singular point that is persistent to variations of the parameters 
may be a symptom that the model is overdetermined in a way that is not obvious to cursory inspection. The singular 
point is defined by the bifurcation problem — the stationary-state equation of the dynamical system — plus equations 
' for the zeros of certain derivatives of the bifurcation problem. This augmented system may have more equations than 
: unknowns because one or more terms incorporating additional parameters are missing. In the language of singularity 
' theory ^ , the codimension (see below for definition) of a persistent, degenerate singular point exceeds the number of 
• • , auxiliary parameters. An idealized model containing this type of point cannot exhibit the qualitative features of a 
more realistic model where perturbational terms unfold the singularity. What is worse is that real-world experiments, 
where perturbations are inevitably present, are likely to exhibit behavior that cannot be predicted by such a model. 
^ ' The two models investigated here describe L- and H-mode dynamics and ELMs in a unified manner, and were 
derived independently by Sugama and Horton ||] (SH) and Lebedev et al [|j (LDGC). Both models describe the 
coupled evolution of state variables related to the pressure gradient, the shear of the poloidal flow, and the level 
of magnetohydrodynamic fluctuations in the edge region of a tokamak. Since we are concerned mainly with the 
stationary states we do not reproduce the dynamical equations, although it should be kept in mind that the stability 
analysis (which is summarized in the bifurcation diagrams) necessarily refers to the dynamics. In this paper we 
show that a canonical analysis of bifurcations innate to these systems as given provides internal evidence that the 
derivations may have neglected important physics. 

The singularity theory analysis essentially consists of three steps. (1) Each model is formulated as the steady-state, 
scalar bifurcation problem g {x, A) = 0, where x is the chosen state variable and A is the chosen control parameter. The 
bifurcation diagrams are found to contain one or more persistent degenerate singularities. (2) We show that the g are 
(locally) strongly equivalent to simple, generic normal forms h. This solves the following recognition problem: what 
conditions must a given g satisfy in order to evince qualitative equivalence to a given normal form hi Concomitantly, 
we obtain two valuable pieces of information: the character of the most degenerate singularity in each model and 
the codimension of this singularity, defined by the minimum number k of independently variable auxiliary parameters 
required to net all possible qualitative behaviors and obtain a universal unfolding. (3) The bifurcation sets are 
perturbed to obtain universal unfoldings of the form G{x,\,ai, . . . ,a]^) — 0, where the k auxiliary or unfolding 
parameters ai, . . . , are non- redundant and all other unfoldings of g may be extracted from G. Singularity theory 
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is concerned with the quahties of steady-state bifurcation problems that determine the dynamics of an associated 
physical system. The key concepts of codimension and qualitative equivalence, together with the universal unfoldings 
and stability considerations, allow us to construct a complete catalogue of the bifurcation behavior. 
The SH model: This may be expressed as the dimensionless bifurcation equation 

g{u,q,da) = [qdaU^'^ - l) {-q + um {u)) , , , 

m{u) = vP{b + au^-P), ^ ' 

obtained by eliminating in the steady state the two other dynamical variables / and k in favor of u oc the potential 
energy of the pressure gradient. The control parameter q is the power input, da is the reciprocal of the anomalous 
diffusivity, and m(u) is the anomalous viscosity. In Sugama and Horton's numerical work da was set to 1, p was 
given values of —3/2 (case A) and —1 (case B), and a and h were given as positive numerical factors. [Note: The 
dynamical equations also contain a parameter c which cancels from Eq. (|l|).) Figure |l| shows the bifurcation diagrams 
for case A and case B. (In all diagrams stable solutions are indicated by solid lines, unstable solutions by dashed lines, 
and branches of limit cycles by dotted lines marking the maximum and minimum amplitude.) It was assumed that 
the transition from the lower stable solution branch (L-mode) to the upper stable branch (H-mode) must occur at 
the singular point A, where the steady-state shear flow kinetic energy / = (m^ — daq) /cu becomes unphysical. The 
transition is discontinuous for case A and continuous for case B. The H-mode branch becomes unstable at a Hopf 
bifurcation to stable limit cycles, identified as ELMs. The SH model thus predicts hysteresis of the L-H transition 
and oscillating and quiescent H-modes, which accords with recent experimental observations llo|-|l^]. However, the 
derivative discontinuity at A is problematic. For case A the transition was described as first-order, but it occurs at 
what appears to be a highly degenerate point. For case B the transition was described as second-order. It should also 
be noted that the singular point A is persistent to variations in da, a, and b. For these reasons we suspect that there 
may not be enough independent parameters in the model. Solution of the recognition problem, step (2), indicates 
that the model may be overdetermined as a bifurcation problem. 

Proposition 1. — Equation with da = dao, P < ~1 is a germ that is strongly equivalent to the normal form 

h{x,\) = -x^ + Xx. (2) 

(The term "germ" is explained as follows: two functions gi{x, A) and 52(2^, A) are equal as germs if they coincide on 
some neighborhood of a fixed point a;o, Aq.) Proof. — We apply the following theorem, adapted from j^]: Theorem. — A 
germ g{x, A) is strongly equivalent to the normal form h{x, A) = ex^ + 6Xx if and only if, at the fixed point {xo, Aq), 

9 = 9x = 9xx = 9\ = 0, gxxx 7^ 0, g\x 7^ (3) 

where e = sgng^xx, 5 — sgng^a;- In Eq. we identify the state variable u = x and the distinguished parameter q = \ 
and evaluate the defining and non-degeneracy conditions (^ at the point A = [uq, go)- We find that g = 9u ~ 9q ~ 0, 
gnu = 4a(-l-|-p) -A{l+p)/da = for da = dao = (l+p) /a{-l+p), and = 12 (1-l-p) /uodao, 9uq = S/uq. 
Equation (||) for the normal form is inferred. It is the prototypic pitchfork a codimension 2 bifurcation which 
requires two auxiliary parameters for an unfolding that contains, to qualitative equivalence, all possible perturbations 
of 5. We see that the defining conditions for the point A yield a system of /our algebraic equations in what is effectively 
only two variables — u and q. To resolve the overdetermination we propose a universal unfolding of Eq. (^. 
Proposition 2. — The bifurcation function 

G{u,q,da,a) ^ g{u,q,da) + a (4) 

is a universal unfolding of the germ (^ for p < — 1. Jt is equivalent to the prototypic universal unfolding of the 
pitchfork G{x, A, a, j3) = —x^ -f Px"^ + \x + a, where da = daO i (3. The proof is not presented here; instead we 
focus on the qualitative consequences. (The physical interpretation of the unfolding parameter a is discussed below.) 
Specifically, Eq. (^ encapsulates the generic behavior of the SH system. The four qualitatively distinct bifurcation 
diagrams are shown in Fig. ^(a)-(d), of which (a) and (b) are physically relevant because a < leads to dynamical 
violation of the condition / > 0. In (a) the L-H and H-L transitions occur at non-degenerate limit points. No marked 
transition to H-mode occurs at all in (b). Now it can be seen why the unperturbed bifurcation set. Fig. |^(c), and the 
partially perturbed bifurcation set. Fig. ^a), cannot predict the results of experiments. The singularity that exists 
in these sets (point A) is not even present when a is nonzero. We also see that changes in the auxiliary parameters 
around the critical values can lead to incomparably different bifurcation behavior. 

What of case B? Proposition 3. — Equation (|I|) with p — —1 is a germ that is strongly equivalent to the normal 
form /i(x. A) = —x^ + A^, a codimension 1 bifurcation known as the transcritical bifurcation. Proof. — We apply the 
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following theorem from Q : Theorem 2. — A germ g{x,X) is strongly equivalent to the normal form h{x, X) = e (x^ ~ ■^^) 
if and only if, at the fixed point {xq, Aq), 

5 = .9. - .9A - 0, .g.. ^ 0, det f '^^^ ) < (5) 

where e = sgng^x- These conditions in Eq. (|l|) yield g = Qu — Oq — 0, guu = —8a, det (Pg = —4{ada — l)^ /v?. 
Equation (Q) in this special, fragile case is a one-parameter universal unfolding, indifferent to the value of da- It yields 
two qualitatively distinct bifurcation diagrams, shown in Fig. ^. Note that the bifurcation structure here excludes 
the possibility of hysteresis. 

The LDGC model: The steady states are summarized in the bifurcation diagram of Fig. ^, where the control 
parameter (f) is the particle flux and p is the pressure gradient. The lower stable branch is identified as L-mode. At A 
the transition to the intermediate stable branch AB, identified as H-mode, is described by Lebedev et al as analogous 
to a second-order phase transition. At B the system moves onto the p = 1 branch in another continuous transition, 
but is said to remain in H-modc. The first Hopf bifurcation initiates a branch of unstable limit cycles and the second 
terminates a branch of stable limit cycles, identified as ELMs. The point C is the intersection of the p — I branch 
and the unstable AC branch. Near B and C the bifurcation equations may be written, respectively, as 

5b(p, (/)) = 7 (0 - dfJ-p'j (p-l) ^p(d- dm'j , and (6) 
gcip, (t>)^l [p^d - 0^ (p - 1) y/ pdm- (7) 

As before, we use the singularity theory analysis to focus on qualitative structure. Using theorem 2 we find that at 
points B and C there is a transcritical bifurcation, which requires the single auxiliary parameter a' for a universal 
unfolding. The two qualitatively distinct bifurcation diagrams are shown in Fig. |^. In (a) a branch of stable limit 
cycles connects the two stable stationary branches. In (b) the branches of stable stationary solutions are unconnected. 
The structure of the limit cycle branch implies that (on a phase plane) a stable orbit is surrounded by an unstable 
orbit. The point A in Fig. ^ clearly is not unfolded by the one-parameter perturbation. Somewhat surprisingly, it is 
the limit point of the branch of the branch CAB, which actually coincides along AB with the continuous branch QAB. 
(This result is detailed elsewhere.) A limit point is its own universal unfolding, i.e., persistent to small perturbations. 

In summary: (1) The SH model in general is a codimension 2 bifurcation problem, containing a pitchfork, that 
requires two unfolding parameters for a universal unfolding and hence complete determination. The critical values 
of the unfolding parameters a and da are respectively and (1 + p) /a (—1 p), p < —1. (2) The LDGC model is a 
codimension 1 bifurcation problem, containing two transcritical bifurcations. A universal unfolding is provided by a 
single auxiliary parameter a' . The two models are therefore structurally dissimilar in general form. However, the fact 
that they describe the same phenomena suggests that the LDGC model may be a partially collapsed codimension 2 
system, and in a forthcoming work we show that this is indeed the case. A fortiori we can also say that second-order 
phase transitions in these systems, if they exist, could only be observed on variation of at least two parameters 
simultaneously. In many bifurcation problems pitchforks occur in the presence of Z2 equivariance in the governing 
equations for the system, that manifests as a physically invariant property. (A function 4>{x) has Z2 symmetry if 
4>{—x) = —(j){x).) A symmetry arises in the dynamical equations for the SH model because the shear of the poloidal 
flow v' is invariant under the transformation v' — > —v'. The unfolding parameter a can therefore be interpreted as a 
symmetry-breaking term, representing an intrinsic energy (or angular momentum) generation rate that occurs even 
in a pressure gradient of zero. The Z2 invariance of the flow shear is not evident in the bifurcation structure of the 
LDGC model, and a' represents a perturbation of the MHD turbulence level. 

Other models for L-H transitions that have multiple solutions include those where the flow shear is due to ion- 
orbit losses on the plasma edge [|l3|,0 or magnetic field ripple induced particle flux in the core [|5|. We feel that 
singularity theory could play an important role in developing and unifying these models and elucidating the physics of 
L-H transitions. There is a reasonable expectation that different models, if they appeal to the same general physics, 
should belong to the same qualitative universality class even though they may differ quantitatively. A wider question 
is whether the dynamics of infinite-dimensional systems can be approximated by low-dimensional systems such as 
these. The practical advantages are obvious, and developments in inertial manifold theory have shown that the 
long-time-scale behavior of infinite-dimensional dissipative systems can occur in a defined finite-dimensional subspace. 
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FIG. 1. Bifurcation diagrams of the original SH model, da = 1, a = 0.05, b = 0.95, c = 5. (a) case A: p = —3/2, (b) case B: 
p = —1. Labels indicate the sign of the shear flow energy, thus a minus sign means that the branch is unphysical. 




FIG. 2. Bifurcation catalogue for the general SH model, p < —1. (a)-(d) The perturbed diagrams, (a) a = 0.01, da = 1, 
(b) a = 0.01, da = 10, (c) a = -0.01, da = 1, (d) a = -0.01, da = 10. (e) The unperturbed diagram, da = 4, a = 0. a = 0.05, 
b = 0.95, c = 5, p = —3/2. (In (b), (d), and (e) the upper Hopf bifurcation is off-scale.) pd: period-doubling bifurcation. 
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FIG. 3. Bifurcation diagrams for the perturbed SH model, case B. p = —1, da — 1, a = 0.05, b — 0.95, c = 5. (a.) a — 0.01, 
(b) a = -0.01. 
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FIG. 4. Bifurcation diagram for the original LDGC model, d — 0.1, dm ~ 0.05, fj, — 0.25, 7 = 5. 
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FIG. 5. The two bifurcation diagrams for the universal unfolding of the LDGC model, (a) a' — 0.01, (b) a' = —0.01. Other 
parameters as for Fig. ^ 
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